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We use first-principles total-energy calculations based on density functional theory to study 
the site occupancy and magnetic properties of Al-substituted M-type strontium hexaferrite 
SrFei 2 -a;AKOi 9 with x = 0.5 and x = 1.0. We find that the non-magnetic AP"*" ions preferen¬ 
tially replace Fe^"*" ions at two of the majority spin sites, 2a and 12fc, eliminating their positive 
contribution to the total magnetization causing the saturation magnetization Mg to be reduced as 
A1 concentration x is increased. Our formation probability analysis further provides the explanation 
for increased magnetic anisotropy field when the fraction of A1 is increased. Although AP"*" ions 
preferentially occupy the 2a sites at a low temperature, the occupation probability of the 12fc site 
increases with the rise of the temperature. At a typical annealing temperature (> 700 °C) AP"*" 
ions are much more likely to occupy the 12k site than the 2a site. Although this causes the mag¬ 
netocrystalline anisotropy Ki to be reduced slightly, the reduction in Ms is much more significant. 

Their combined effect causes the anisotropy field Ha to increase as the fraction of A1 is increased, 
consistent with recent experimental measurements. 


I. INTRODUCTION 

Strontium hexaferrite, SrFei20i9 (SFO) is one of the 
most commonly used materials for permanent magnets, 
magnetic recording and data storage, and components in 
electrical devices operating at microwave/GHz frequen¬ 
cies, due to its high Gurie temperature, large saturation 
magnetization, excellent chemical stability and low man¬ 
ufacturing cost m- However, in comparison with Nd- 
Fe-B and magnet, the coercivity of the SFO is low and 
presents a significant limitation in its application. There¬ 
fore, enhancement of the coercivity is an important re¬ 
search topic for the strontium hexaferrite. 

In order to tailor the magnetic properties such as mag¬ 
netization and coercivity, various cation substitutions in 
the M-type hexaferrites have been investigated. For ex¬ 
ample, the substitution of La [3 [7], Sm i, Pr (2 and 
Nd [To] in the SFO increased coercivity moderately while 
the substitution of Zn-Nb nn, Zn-Sn [TTHTT] and Sn-Mg 
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m decreased coercivity. However, the coercivity of the 
M-type hexaferrites is not increased significantly by these 
cation substitutions, and is still much smaller than that 
of Nd-Fe-B magnet [T5] . 

A1 substitution in the M-type hexaferrite has been 
more effective in enhancing coercivity [I6l - l20] . Par¬ 
ticularly, Wang et al synthesized Al-doped SFO 
SrFei 2 -xAlxOi 9 (SFAO) with Al content of a; = 0 — 4 
using glycinnitrate method and subsequent annealing in 
a temperature over 700 °C obtaining the largest coer¬ 
civity of 17.570 kOe, which is much larger than that of 
SFO (5.356 kOe) and exceeds even the coercivity of the 
NdaFeirB (15.072 kOe) p. Wang and co-workers also 
observed that the coercivity of the SFAO increases with 
increasing Al concentration at a fixed annealing temper¬ 
ature. These results call for a systematic understand¬ 
ing, from first principles, of why certain combinations of 
dopants lead to particular results. This theoretical un¬ 
derstanding will be essential in systematically tailoring 
the properties of SFO. 

There have been several previous first-principles inves¬ 
tigations of SFO. Fang et al investigated the electronic 
structure of SFO using density-functional theory (DFT) 
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m- Park et al have calculated the exchange interaction 
of SFO from the differences of the total energy of dif¬ 
ferent collinear spin configurations [12]. In spite of the 
importance of substituted SFO, only a few theoretical 
investigations have been done. Magnetism in La substi¬ 
tuted SFO has been studied using DFT |23l |2l|. The site 
occupancy and magnetic properties of Zn-Sn substituted 
SFO has been investigated El- 
In this work we use first-principles total-energy calcu¬ 
lations to study the site occupation and magnetic prop¬ 
erties of Al substituted M-type strontium hexaferrite 
SrFei 2 -xAla;Oi 9 with x = 0.5 and x = I.O. Based on 
DFT calculations, we determine the the structure of var¬ 
ious configurations of SFAO with different Al concentra¬ 
tions and compute the occupation probabilities for dif¬ 
ferent substitution sites at elevated temperatures. We 
show that our model predicts an decrease of saturation 
magnetization as well as a decrease in magnetocrystalline 
anisotropy ATi, and the increase of the anisotropy field Ha 
as the fraction of Al is increased, consistent with recent 
experimental measurements. 



II. METHODS 

SFO has a hexagonal magnetoplumbite crystal struc¬ 
ture that belongs to P63/totoc space group. Fig.[^shows 
a unit cell of SFO used in the present work that con¬ 
tains 64 atoms of two formula units. Magnetism in SFO 
arises from Fe^"*" ions occupying five crystallographically 
inequivalent sites in the unit cell, three octahedral sites 
(2a, 12fc, and 4 / 2 ), one tetrahedral site (4/i), and one 
trigonal bipyramidal site {2b) as represented by the poly- 
hedra in Fig. Ea). SFO is also a ferrimagnetic material 
that has 16 Fe^ ions with spins in the majority direction 
(2a, 26, and 12k sites) and 8 Fe^'*' ions with spins in the 
minority direction (4/i and 4/2 sites) as indicated by the 
arrows in Fig. [^b). 

Total energies and forces were calculated using DFT 
with projector augmented wave (PAW) potentials as im¬ 
plemented in VASP jlSl US] . All calculations were spin 
polarized according to the ferrimagnetic ordering of Fe 
spins as first proposed by Gorter iniiiT]. A plane-wave 
energy cutoff of 520 eV was used both for pure SFO and 
Al-substituted SFO. Reciprocal space was sampled with 
a 7 X 7 X 1 Monkhorst-Pack mesh [IS] with a Fermi- 
level smearing of 0.2 eV applied through the Methfessel- 
Paxton method [29]. We performed relaxation of the elec¬ 
tronic degrees of freedom until the change in free energy 
and the band structure energy was less than 10“^ eV. We 
performed geometric optimizations to relax the positions 
of ions and cell shape until the change in total energy 
between two ionic step was less than 10“"^ eV. Electron 
exchange and correlation was treated with the general¬ 
ized gradient approximation (GGA) as parameterized by 
the Perdew-Burke-Ernzerhof (PBE) scheme [30]. To im¬ 
prove the description of localized Fe 3d electrons, we em¬ 
ployed the GGA-I-U method in the simplified rotationally 


FIG. 1. (color online) (a) A unit cell of SFO containing 
two formula units. Two large gold spheres are Sr atoms and 
small gray spheres are O atoms. Colored spheres enclosed 
by polyhedra formed by O atoms represent Fe®^ ions in dif¬ 
ferent inequivalent sites: 2a (blue), 26 (cyan), 12fc (purple), 
4/1 (green), and 4/2 (red), (b) A schematic diagram of the 
lowest-energy spin configuration of Fe®”*" ions of SFO. The 
arrows represent the local magnetic moment at each atomic 
site. (For interpretation of the references to color in this fig¬ 
ure caption, the reader is referred to the online version of this 
paper.) 

invariant approach described by Dudarev et al |31j . This 
method requires an effective U value {Ues) equal to the 
difference between the Hubbard parameter U and the ex¬ 
change parameter J. We chose UeS equal to 3.7 eV for 
Fe based on the previous result [H] . 

III. RESULTS AND DISCUSSION 

The substitution of Fe^"*" ions by Al^+ ions consider¬ 
ably affects the unit cell parameters. We have calculated 
the lattice parameters of pure and Al-substituted SFO by 
relaxing ionic positions as well as the volume and shape 
of the unit cell. In all cases the final unit cell was found 
to remain hexagonal. In the case of pure SFO, the lat¬ 
tice parameters a and c were found to be 5.93 A and 
23.21 A in good agreement with the experimental values 
of a = 5.88 A and c = 23.04 A, respectively [TU 15^ : 
the deviation between the experimental and the theo¬ 
retical values is less than 1%. In the case of a; = 0.5 
in SrFei 2 -a;Al 2 ,Oi 9 the lattice parameters a and c were 
calculated to be 5.92 A and 23.16 A respectively, while 
the volnme of the unit cell was reduced by 0.61%. For 
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FIG. 2. (color online) Comparison of calculated and experi¬ 
mental (Ref. [19]) volume of the unit cell of SrFei 2 -i,Ala;Oi 9 
as a function of the fraction of A1 x. 


X = 1.0, a = 5.91 A and c = 23.04 A were found, and 
reduction in the unit cell volume was 2.51%. Fig.j^shows 
that the reduction of unit cell volume predicted by our 
DFT calculation is consistent with the experimental re¬ 
sults [Hill]. 

We investigated the site preference of A1 substituting 
Fe in SrFei 2 -a;Ala;Oi 9 for (i) x = 0.5 and, (ii) x = 1.0. 
The X = 0.5 case corresponds to the condition where 
one A1 atom is substituted in the unit cell, while two A1 
atoms were substituted in the case of cc = 1.0 as shown 
Fig.i To determine the site preference of the substituted 
A1 atoms, the substitution energy of configuration i was 
calculated using the following expression: 

F;sub(») = £;(SFAO(^)) - F;(SF0) - ^ n„e(a) (1) 


where A(SFAO(f)) is the total energy per unit cell at 
0 K for SFAO in configuration i while A(SFO) is the 
total energy per unit cell at 0 K for SFO. e{a) is the 
total energy per atom for element a {a = Al, Fe) at 0 K 
in its most stable crystal structure, ria is the number 
of atoms of type a added: if two atoms are added then 
Ua = +2 while = — 1 when one atom is removed. 
The configuration with the lowest Eguh is concluded to 
be the ground state configuration, and the corresponding 
substitution site is the preferred site of Al atoms at 0 K. 

To understand the site preference of the substituted 
Al3+ ions at higher temperatures, we compute the for¬ 
mation probability of configuration i using the Maxwell- 
Boltzmann statistical distribution [33] : 

p ^ g»exp(-AG,/fcBT) 

Ej 9j exp(-ACj/fesT) 



FIG. 3. (color online) The structures of SrFei 2 - 2 ;Ala;Oi 9 
with spins oriented in the easy axis (001): (a) configuration 
[2a] for X = 0.5 and (b) configuration [2a, 12fc].l for x = 1.0. 
Al atoms are labeled and other atoms are colored as in Fig.[^ 


where gi is the multiplicity of configuration i (number of 
equivalent configurations) and 


AG, = AE,uh{i) + PAV, - TAS, (3) 

is the change of the free energy of configuration i rela¬ 
tive to that of the ground state configuration; AAsub(*)j 
AVi, and ASi are the substitution energy change, volume 
change, entropy change for configuration 1; P and fee are 
the pressure and Boltzmann constant. 

For the x = 0.5 concentration, one Al atom is substi¬ 
tuted at one of the 24 Fe sites of the unit cell as shown in 
Fig. j^a). The application of crystallographic symmetry 
operations shows that many of these Fe sites are equiv¬ 
alent and leaves only five inequivalent structures. We 
label these inequivalent configurations using the crystal¬ 
lographic name of the Fe site: [2a], [26], [4/i], [ 4 / 2 ], and 
[12fc]. These structures were created by substituting one 
Al atom to the respective Fe site of a SFO unit cell and 
performing full optimization of the unit cell shape and 
volume, and ionic positions. 

Table HI lists the results of our calculation for all five in¬ 
equivalent configurations in the order of increasing sub¬ 
stitution energy. The lowest Esuh is found for configu¬ 
ration [2a] shown in Fig. [^a). We can conclude that 
at 0 K the most preferred site for the substituted Al 
atom is the 2a site. We used Eq. ([^ to compute the 
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TABLE I. Five inequivalent configurations of 
SrFei 2 -a;ALOi 9 with x = 0.5. g is the multiplicity of 
the configuration. i?sub is the substitution energy of the 
SFAO. The total magnetic moment (mtot) and its change 
with respect to SFO (Amtot) are also given. All values are 
for a double formula unit cell containing 64 atoms. Energies 
are in eV while magnetic moments are in ^b. 


Config 

9 

Esuh 

mtot 

Amtot 

[2a] 

2 

-6.04 

35 

-5 

[12fc] 

12 

-6.00 

35 

-5 

[ 4 / 2 ] 

4 

-5.63 

45 

-f5 

[26] 

2 

-5.60 

35 

-5 

[4/ij 

4 

-5.57 

45 

-f5 



FIG. 4. (color online) Temperature dependence of the forma¬ 
tion probability of different configurations of SrFei 2 -iAla;Oi 9 
with X = 0.5. The configurations with negligible probability 
are not shown. 


probability to form each configuration as a function of 
temperature. Since the volume change among different 
configurations is very small (less than 0.1 we can 
safely regard PAV term to be negligible (in the order 
of 10“^ eV at the standard pressure of 1 atm) com¬ 
pared to the Aiiisub(*) term in Eq. (§. The entropy 
change AS' has a configurational part, ASc, and a vi¬ 
brational part, ASvib |34j . For binary substitutional al¬ 
loys such as the present system, ASvib is around 0.1- 
0.2 AiB/atom, and ASc is 0.1732 /ce/atom [33]. Therefore, 
we set AS = 0.3732 fee/atom. 


Fig. HI displays the temperature dependence of the 
formation probability of different configurations of 
SrFei 2 -a;Ala;Oi 9 with x = 0.5. The doped Al^+ ions 
mainly replace Fe^+ ions from the 2a and the 12fc sites. 


The formation probabilities of 
negligible and not shown in Fig. 
the doped Al^+ ion replaces Fe^ 


26], [4/i] and [ 4 / 2 ] are 
^ The probability that 
ion from the 2a site is 


maximum at 0 K and it falls as temperature increases. 


TABLE 11. Ten lowest energy inequivalent configurations of 
SrFei 2 -a;AbOi 9 with x = 1.0. g is the multiplicity of the con¬ 
figuration. Eaub is the substitution energy per A1 atom. The 
total magnetic moment (mtot) and its change with respect to 
SFO (Amtot) are also given. All values are for a double for¬ 
mula unit cell containing 64 atoms. Energies are in eV while 
moments are in fiB- 


Config 

9 

Esuh 

mtot 

Amtot 

[2a, 2a] 

1 

-6.056 

30 

-10 

[2a, 12fc].l 

12 

-6.054 

30 

-10 

[2a, 12fc].2 

12 

-6.041 

30 

-10 

[12fc, 12fc].l 

6 

-6.025 

30 

-10 

[12fe, 12fc].2 

12 

-6.025 

30 

-10 

[12fe, 12fc].3 

12 

-6.027 

30 

-10 

[12fc, 12fc].4 

12 

-6.025 

30 

-10 

[12fc, 12fc].5 

6 

-6.023 

30 

-10 

[12fc, 12fc].6 

6 

-6.017 

30 

-10 

[12fe, 12fc].7 

12 

-6.014 

30 

-10 


while the occupancy of Al^+ at the 12fc site rises with 
temperature. The two curves cross at T ~ 220 K. At 
a typical annealing temperature of 1000 K for SFAO [1] 
the site occupation probability of the site 2a and 12fc is 
0.196 and 0.798, respectively. Thus, during the anneal¬ 
ing process of the synthesis of the SFAO the doped Al^+ 
ions are more likely to replace Fe^’*' ions from the 12fc site 
than the 2a site despite of higher substitution energy. 

For the a; = 1.0 concentration, two A1 atoms are substi¬ 
tuted at two of the 24 Fe sites of the unit cell as shown in 
Fig.§b). These Fe sites have more than one equivalent 
site. Substitution of A1 atoms breaks the symmetry of the 
equivalent sites of pure SFO. Out of all C(24, 2) = 276 
possible structures, many of the structures are crystal- 
lographically equivalent. On applying crystallographic 
symmetry operations, the number of inequivalent struc¬ 
tures reduces to 40. We label these inequivalent configu¬ 
rations using the convention of [(site for the first Al),(site 
for the 2nd Al)].(unique index). For example, when two 
A1 atoms are substituted at the 2a and 12A: sites, there 
are 2 inequivalent configurations, which are labeled as 
[2a, 12A:].l and [2a, 12A:].2. These structures are fully op¬ 
timized and their substitution energies are calculated us¬ 
ing Eq. 0 - When there are more than one inequivalent 
configuration, we assign the unique index in the order of 
increasing Esuh- 

Table [^ lists the ten lowest energy configurations of 
SrFei 2 -a:Ala;Oi 9 with x = 1.0. The configuration [2a, 2a] 
where two Al^+ ions replace Fe^"*" ions from two 2a sites 
has the lowest Fgub, and it is the most energetically fa¬ 
vorable configuration at 0 K. To investigate the site oc¬ 
cupation at nonzero temperatures we compute the for¬ 
mation probability of each configuration using Eq. (|^. 
Similar to the previous case the volume change among 
different configurations is very small (less than 0.7 
and we can safely ignore the PAV term. The entropy 
term is calculated in the same way as the x = 0.5 case. 
Fig.0 shows the variation of the formation probability 
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FIG. 5. (color online) Temperature dependence of the forma¬ 
tion probability of different configurations of SrFei 2 - 2 ;Ala;Oi 9 
with X = 1.0. For clarity only the configurations with signifi¬ 
cant formation probability are labeled. 


of different configurations with temperature. We note 
that due to low multiplicity of the configuration [2o, 2a], 
its formation probability falls rapidly as temperature in¬ 
creases. On the other hand, the formation probability of 
the configuration [2a, 12fc] (sum of the formation proba¬ 
bilities for all [2a, 12k].n configurations) increases steeply 
and reaches a maximum value at 50 K and then falls with 
temperature. Fig.j^shows that the formation probability 
of the [2o, 12fc] conhguration becomes larger that that of 
[2a, 2a] beyond T ~ 10 K, which is a much lower transi¬ 
tion temperature than in the x = 0.5 case. 

We can calculate the occupation probability of A1 at 
nonzero temperatures for a given site by adding all for¬ 
mation probabilities of the configurations where at least 
one Al^+ ion is substituted in that site. At the annealing 
temperature of 1000 K, the occupation probability of A1 


for 12fc site is 79.8% for x = 0.5 as given in Table IV 


The same probability is increased to 97.7% for x = 1.0 
as calculated by adding the Piooo’s for all configurations 
that contain the 12fc site. This means that the fraction 
of AP+ ions occupying the 12fc site increases when the 
fraction of A1 is increased from x = 0.5 to a:: = 1.0. This 
conclusion is in agreement with the previously reported 
measurements mill 135]. 


In Table |III| we compare the contribution of differ¬ 
ent sublattices to the total magnetic moment in Al- 
substituted SFO. To see the effect of Al^+ ions in differ¬ 
ent substitution sites, we split the entries of sublattices 
containing these ions (2a and 12fc). As expected, Al^+ 
ions carry negligible magnetic moment regardless of their 
substitution sites. Consequently, when they replace Fe^+ 
ions in the minority spin sites (4/i and 4 / 2 ), they elimi¬ 
nate a negative contribution and hence increase the total 
magnetic moment. On the other hand, when they re¬ 


place Fe^’*' ions in the majority spin sites (12fc, 2a, and 
26), they eliminate a positive contribution and hence re¬ 
duce the total magnetic moment. For the x = 0.5 case, 
the most probable sites are 12fc and 2a (majority sites) 
and the net magnetic moment of the unit cell is reduced 
by 5 For the configuration [2a, I2A:].l of the x = 1.0 
case, two A1 atoms are substituted in the 2o and 12A: 
sites, there is a reduction of 10 /tb in the total magnetic 
moment per unit cell. 

Magnetic Anisotropy determines the capacity of a 
magnet to withstand external magnetic and electric 
fields. This property is of considerable practical in¬ 
terest, because anisotropy is exploited in the design of 
the most magnetic materials of commercial importance. 
The magnetocrystalline anisotropy energy (MAE) is one 
of the main factors that determine the total magnetic 
anisotropy of the material. To investigate the effect of 
A1 substitution on the magnetic anisotropy of SFO, we 
computed the MAE and the magnetic anisotropy con¬ 
stant of SrFei 2 -a;Ala;Oi 9 for x = 0, 0.5 and 1. The MAE, 
in the present case, is defined as the difference between 
the two total energies where electron spins are aligned 
along two different directions |36j : 


Amae — A(ioo) — A(ooi) 


(4) 


where A(ioo) is the total energy with spin quantization 
axis in the magnetically hard plane and A(ooi) is the to¬ 
tal energy with spin quantization axis in the magneti¬ 
cally easy axis. Using the MAE, the uniaxial magnetic 
anisotropy constant Ki can be computed [33I3H1 


Ki 


Emae 
V sin^ 0 


(5) 


where V is the equilibrium volume of the unit cell and 
0 is the angle between the two spin quantization axis 
orientations (90° in the present case). The anisotropy 
field Ha can be expressed as [33] 


Ha 


2Ki 


( 6 ) 


where Ki is a magnetocrystalline anisotropy constant 
and Ms is saturation magnetization. 

The results for the MAE, the magnetocrystalline 
anisotropy constant Ki, and anisotropy field Ha for 
SFAO with different A1 concentration are presented in 
Table |IV[ To compare with experimental results, we also 
compute the weighted average of Ki and Ha using the 
formation probability Piooo at a typical annealing tem¬ 
perature of 1000 K [T]. We note that SFAO considered in 
the present work loses most of its magnetic properties at 
typical annealing temperatures (1000 K or higher) that 
are near or above its Curie temperature. The magnetic 
properties listed in Table W refer to their ground state 
properties at the temperature T = 0. We use the for¬ 
mation probability at 1000 K to compute the weighted 
averages as the crystalline conhgurations of SFAO will be 
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TABLE III. Contribution of atoms in each sublattice to the total magnetic moment of Al-substituted SFO structnres [12A:], 
[2a], and [2a, 12A:].l compared with pure SFO. All magnetic moments are in /tb. Am is measnred relative to the values for 
the pure SFO. Note that the total magnetic moment of the unit cell (mtot) is slightly different than the snm of local magnetic 
moments due to the contribution from the interstitial region. 


site 

SFO 

[12fc] 

[2a] 

[2a, 12k\.l 

atoms 

m 

atoms 

m 

Am 

atoms 

m 

Am 

atoms 

m 

Am 

2d 

2 Sr 

-0.006 

2 Sr 

-0.006 

0.000 

2 Sr 

-0.006 

0.000 

2 Sr 

-0.006 

0.000 

2a 

1 Fe 

4.156 

1 Fe 

4.155 

-0.001 

1 A1 

-0.010 

-4.166 

1 A1 

-0.010 

-4.166 

1 Fe 

4.156 

1 Fe 

4.156 

0.000 

1 Fe 

4.156 

0.000 

1 Fe 

4.156 

0.000 

2b 

2 Fe 

8.098 

2 Fe 

8.086 

-0.012 

2 Fe 

8.100 

0.001 

2 Fe 

8.090 

-0.008 

4/i 

4 Fe 

-16.152 

4 Fe 

-16.189 

-0.037 

4 Fe 

-16.268 

-0.116 

4 Fe 

-16.304 

-0.152 

4/2 

4 Fe 

-16.384 

4 Fe 

-16.420 

-0.036 

4 Fe 

-16.382 

0.002 

4 Fe 

-16.418 

-0.034 

12fc 

1 Fe 

4.172 

1 A1 

0.000 

-4.172 

1 Fe 

4.170 

-0.002 

1 A1 

-0.001 

-4.173 

11 Fe 

45.884 

11 Fe 

45.861 

-0.023 

11 Fe 

45.872 

-0.012 

11 Fe 

45.846 

-0.038 

4e 

4 0 

1.416 

4 0 

1.304 

-0.112 

4 0 

1.424 

0.008 

4 0 

1.316 

-0.100 

4 / 

4 0 

0.360 

4 0 

0.281 

-0.079 

4 0 

0.310 

-0.050 

4 0 

0.230 

-0.129 

6h 

6 0 

0.124 

6 0 

0.115 

0.009 

6 0 

0.134 

0.010 

6 0 

0.117 

-0.007 

I2k 

12 0 

1.016 

12 0 

0.877 

-0.129 

12 0 

0.548 

-0.468 

12 0 

0.404 

-0.612 

12k 

12 0 

2.140 

12 0 

1.895 

-0.245 

12 0 

2.088 

-0.052 

12 0 

1.839 

-0.301 



38.980 


34.114 

-4.837 


34.136 

-4.845 


29.259 

-9.720 

mtot 


40 


35 

-5 


35 

-5 


30 

-10 


TABLE IV. The satnration magnetization (Ms), magnetocrystalline anisotropy energy (MAE), magnetocrystalline anisotropy 
constant (Ki) and anisotropy field (Ha) for SFO and SFAO. x is the A1 fraction in SrFei2-2,AL0ig and V is the volume of the 
unit cell in A^. Piooo is the formation probability at 1000 K. The averaged quantities are weighted by Piooo- Ms is in emu/g, 
MAE in meV, Ha in kOe, and Ki in kj-m”®. 


X 

Config 

Ms 

MAE 

V 

Ki 

Ha 

Piooo 

{Ms) 

{Ki) 

{Ha) 

0.0 

SFO 

110.19 

0.85 

707.29 

193 

7.35 

1.000 

110.19 

193 

7.35 


[2a] 

96.41 

0.95 

703.29 

216 

9.38 

0.196 





[12k] 

96.41 

0.80 

703.19 

182 

7.90 

0.798 




0.5 

[2b] 

96.41 

0.67 

702.82 

152 

6.62 

0.003 

96.49 

189 

8.18 


[ 4 / 1 ] 

123.96 

0.86 

704.22 

196 

6.61 

0.001 





[ 4 / 2 ] 

123.96 

0.83 

702.58 

189 

6.38 

0.001 





[2a, 2a] 

82.64 

0.99 

698.94 

227 

11.41 

0.019 





[2a, 12fc] 

82.64 

0.88 

699.08 

202 

10.13 

0.379 





[12fc, 12fc] 

82.64 

0.75 

698.66 

172 

8.64 

0.585 





[12fc,4/2[ 

110.19 

0.78 

690.64 

181 

6.74 

0.007 





[12fc,4/i] 

110.19 

0.80 

700.29 

183 

6.92 

0.004 





[12fc,26] 

82.64 

0.62 

698.98 

142 

7.14 

0.002 





[ 4 / 2 , 4 / 2 ] 

137.74 

0.80 

697.96 

184 

5.53 

0.000 




1.0 

[ 4 / 2 , 4 / 1 ] 

137.74 

0.83 

699.62 

191 

5.74 

0.000 

83.03 

184 

9.23 


[4/2,26] 

110.19 

0.60 

698.82 

138 

5.19 

0.000 





[ 4 / 2 , 2a] 

110.19 

0.90 

698.53 

206 

7.78 

0.002 





[ 4 / 1 , 4 / 1 ] 

137.74 

0.86 

701.38 

196 

5.95 

0.000 





[4/1,26] 

110.19 

0.65 

699.95 

149 

5.62 

0.000 





[ 4 / 1 ,2a] 

110.19 

0.91 

700.11 

208 

7.87 

0.001 





[26, 26] 

82.64 

0.45 

698.86 

103 

5.19 

0.000 





[26, 2a] 

82.64 

0.74 

698.82 

170 

8.53 

0.001 





distributed according to this value during the annealing 
process. 

Table HY] shows that Mg decreases as the concentration 
of A1 X is increased from 0 to 0.5 to 1.0, consistent with 
the previous experimental results m uni Eoi E]. Our 
calculation also shows that Ki decreases as the concen¬ 
tration of A1 X is increased from 0 to 0.5 to 1.0. At a 
low temperature A1 atoms prefer to occupy the 2a sites, 


which would have increased Ki (see Ki values for [2a] 
and [2a, 2a] in Table [TV] ) . However, the formation prob¬ 
ability of the configurations involving 12A: site (such as 
[12fc], [2a, 12fc] and [12A:, 12fc]) increases significantly as 
the temperature rises due to the entropy contribution of 
the free energy. At the annealing temperature Al^+ ions 
are much more likely to occupy the 12fc site than the 2a 
site. This causes the magnetocrystalline anisotropy con- 
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stant Ki of Al-substituted SFO to be reduced with the 
increase of A1 fraction x, consistent with the experimen¬ 
tal measurement reported by Albanes [ID]. Despite of 
this, Ms is reduced more significantly than Ki and this 
causes the anisotropy field Ha in Eq. (|^ to increase as 
the concentration of A1 x is increased from 0 to 0.5 to 
1.0 as shown in Table CYl This result is consistent with 
several experimental results Dill]- 

IV. CONCLUSIONS 

Using the first-principles total energy calculations 
based on density functional theory, we obtained the 
ground state structures and associated formation prob¬ 
abilities at finite temperatures for Al-substituted SFO, 
SrFei 2 -a;AUOi 9 with x = 0.5 and 1.0. The structures de¬ 
rived from our calculations show that the total magnetic 
moment of the SFO unit cell is reduced as the fraction 
of A1 atoms increases. This reduction of magnetization 
is explained by the fact that the non-magnetic A1 atoms 
prefer to replace Fe^"*" ions at two of the majority spin 
sites, 2a and 12fc, eliminating their positive contribution 


to the total magnetization. Our model also explains the 
increase of the observed anisotropy field when the frac¬ 
tion of A1 in SFO is increased. At the annealing temper¬ 
ature AU+ ions are much more likely to occupy the \2k 
site than the 2a site. Although this causes the magne¬ 
tocrystalline anisotropy to decrease slightly, the reduc¬ 
tion in the saturation magnetization is larger and their 
combined effect causes the magnetic anisotropy field of 
Al-substituted SFO to be reduced with increase of A1 
fraction x. Our results are consistent with the available 
experimental measurement on Al-substituted SFO. 
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